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ABSTRACT 


We present a new class of solutions for the gas flows in elliptical galaxies 
containing massive central black holes (BH). Modihed King model galaxies 
are assumed. Two source terms operate: mass loss from evolving stars, and 
a secularly declining heating by supernovae (SNIa). Relevant atomic physical 
processes are modeled in detail. Like the previous models investigated by Ciotti 
et al. (1991, CDPR), these new models hrst evolve through three consecutive 
evolutionary stages: wind, outflow, and inflow. At this point the presence of the 
BH alters dramatically the subsequent evolution, because the energy emitted 
by the BH can heat the surrounding gas to above virial temperatures, causing 
the formation of a hot expanding central bubble. Short and strong nuclear 
bursts of radiation (Lbh) are followed by longer periods during which the X-ray 
galaxy emission comes from the coronal gas (Lx). The range and approximate 
distribution spanned by Lx are found to be in accordance with observations 
of X-ray early type galaxies. Moreover, although high accretion rates occur 
during bursting phases when the central BH has a luminosity characteristic of 
quasars, the total mass accreted is very small when compared to that predicted 
by stationary cooling-flow solutions and computed masses are in accord with 
putative BH nuclear masses. In the bursting phases the X-ray gas luminosity 
is low and the surface brightness prohle is very low compared to pre-burst 
or to cooling flow models. We propose that these new models, while solving 
some long-standing problems of the cooling flow scenario, can provide a unihed 
description of QSO-like objects and X-ray emitting elliptical galaxies, these 
being the same objects observed at two different evolutionary phases. 

Subject headings: galaxies: active — cooling flows — evolution — X-rays 
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1. Introduction 


As first revealed by Einstein observations, normal elliptical galaxies, both isolated or 
in groups and clusters, can be powerful X-ray sources with 0.5 - 4.5 KeV luminosities Lx 
ranging from ~ 10^® to ~ 10^^ erg s“^. This emission is associated with hot gaseous halos 
within the galaxies, containing Mgas = 10® — IO^^Mq (see [Fabbiano 1989|) . 

In order to explain this observational hnding, a certain class of solutions designated 
cooling flow models have been proposed and extensively investigated (e.g., Fabian, Nulsen, 
& Canizares 1984; Sarazin & White 1987,1988; Vedder, Trester, & Canizares 1988). While 
these models have many attractive features, they are far from giving a totally satisfactory 
account of the X-ray properties of all elliptical galaxies, as most observed systems are 
much fainter in the X-rays than the models predict and have different radial prohles than 
expected. Moreover, the cooling flow models do not solve the question of where the cool gas 
is deposited: over a Hubble time an amount of material comparable to the mass of stars in 
the galactic core flows into the nucleus, but the expected distortions of the central optical 
surface brightness and velocity dispersion are not observed. 


One possible solution to part of the previous set of problems was proposed by D’Ercole 
et ah (1989) and CDPR, who showed that the heating from SNIa could be effective in 
maintaining low luminosity galaxies in a wind phase over an Hubble time (and so preventing 
the gas from accumulating in the centre). But the most massive galaxies ultimately 
experience a central cooling catastrophe, leading to a situation similar to a cooling flow. 
Clearly a component of the explanation is missing and is possibly related to the fact (e.g., 
|Rees 198'4D that many (perhaps most) early-type galaxies show a nuclear activity, and, 
according to the standard interpretation of the AGN phenomenon, a massive BH is at its 
origin. So it is natural to investigate the accretion of a galactic gas inflow onto galaxies 
within which lurk massive central BHs (Mbh ~ 1O®M0). 
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Binney & Tabor (1995, BT) explored this problem with the aid of spherically symmetric 
numerical simulations, assuming an homogeneous release of energy in the inner kpc of their 
galaxy models during the inflow phases. Moreover, BT assumed that all the accretion 
luminosity was available for the gas heating, due to the interaction between a nuclear jet 
and the surrounding ISM. 


In the present paper we explore, by numerical integration of the fully non-stationary 
equations of hydrodynamics, the modihcations on the results of CDPR, assuming the 
presence of a massive BH the galaxy centre with detailed allowance for the effects on the 
flow of the radiation emitted by the central BH. As will be shown, the gas over the body 
of the galaxy is (as noted by BT) really optically thin, but nevertheless the effect of energy 
exchange between the nuclear radiation and the gas flow is dramatic. This effect was 
already known and extensively studied for accreting compact objects ( Pstriker et al. 1976|; 
Cowie, Ostriker, & Stark 1978). In a successive paper (|Ciotti fc Pstriker 1997| , Paper II), a 
quantitative analysis of various aspects of the scenario summarized in this Letter, together 
with an exaustive description of the input physics and its modelization, will be given. 


2. Results 

All the results shown here refer to a model whose parameters are hxed following 
the line of CDPR. The stellar density prohle is a King (1972) distribution, with total 
blue luminosity Lb = 5 x lO^^L©, central velocity dispersion Uo* = 280 km s“^, and core 
radius Uc* = 350 pc. The dark-matter halo is described by a quasi-isothermal density 
distribution, with Mh/M* = 7.8 and r^h/Tc* = 4.2. The SNIa rate is the same as that in 
the King Reference Model of CDPR. The bolometric luminosity emitted by accretion onto 
the BH is Lbh = cc^ATbe, where c is the light velocity, and e is the accretion efficiency, 
with 10“^^e^l0“^. The spectral distribution of Lbh near the BH is assumed to be 
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-^bh(^) cc where /iz/b=lMeV. In the present model e = 0.1. 

The spherically symmetric hydrodynamical equations are integrated numerically using 
the Eulerian up-wind scheme with time splitting and artihcal viscosity as used in CDPR. 
In the energy equation the contribution of Compton heating (and cooling) of the gas due 
to Lbh and to the recycling of the bremsstrahlung radiation produced by the gas heated 
by the BH activity are included. We allow also for the effect of the photoionization in 
both cooling and heating of cold gas as well as momentum exchange between photons and 
electrons. At each radius the radiation held is integrated, considering the gas differential 
absorption on Lbh(^) and using for the electrons the Klein-Nishina cross-section. In this 
way the absorbed fraction of Lbh is computed self-consistently. 

In Fig. lb (solid line) the temporal evolution of the coronal X-ray luminosity Lx of 
the gas in the 0.5-4.5 KeV band is shown over an Hubble time. The evolution up to the 
so-called cooling catastrophe (t ~ 9.4 Gyr) is analogous to that described in CDPR, but 
after this time the Compton heating instability completely alters the how evolution and 
its properties. At the cooling catastrophe negative infall velocities appear near the galaxy 
center, with Mbh ~ 60 and this accretion produces a strong, energetic feed-back 

producing a very high Lbh (Fig. la). The gas in the central regions of the galaxy is 
strongly heated to temperatures comparable with the Compton temperature associated 
with Lbh(^) (— 10® K), and starts to expand, decreasing its density by more than two 
orders of magnitude, driving a shock wave outwards and producing a hot bubble of a few 
hundred parsecs in diameter. The net ehect is, observationally, a large reduction of Lx 
(Fig. lb), and, hydrodynamically, the interruption of the galactic inhow and the consequent 
shut-oh of Lbh- Then the radiative losses increase again, and, after a period of the order 
of the hot gas cooling time, the cycle repeats. In the model described here this time is 
of the order of ~ 1 Gyr. In the case of very high accretion the shock wave can reach the 
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galaxy edge, and expel gas from the galaxy. At higher time resolution each burst shows a 
very complex structure, that will be discussed in detail in Paper II; the temporal blow-up 
of the hrst burst shown in Fig. 1 is plotted in Fig. 2, showing QSO-like luminosities. An 
important characteristic of all computed models - of which a single representative is here 
discussed - is that the fraction of Lbh effectively absorbed by the gaseous halo is in the 
range 10“^ — 10“^, but the gas flows are found to he invariably unstable due to Compton 
heating for all the explored efficiencies: in presence of a massive BH at the center of 
elliptical galaxies the possibility of a stationary cooling flow seems to be very remote. 

In Fig. 3 we show the distribution of Lx from 9 Gyr to 15 Gyr. The dashed histogram 
shows the model distribution of Lx given in Fig. 1, and the solid line shows data for 
(non-boxy) early-type galaxies taken from Fig.l of GDPR. Finally the dotted histogram is 
the distribution of Lx for the same model, with the cooling flow assumption of e = 0. We 
see that the model with e = 0.1 has a distribution, over time, of Lx surprisingly similar 
to that of observed galaxies, but the cooling flow model (as is well known) produces far 
too much radiation. In Fig. 4 the X-ray surface brightness prohle (Sx) of the presented 
model is shown at two different epochs, before and during bursts (vertical arrows in Fig. 2). 
Also shown is the cooling flow prohle for the same galaxy at t = 15 Gyr. Note how Sx is 
characterized by a well dehned core before a burst, alleviating the problem of the too cuspy 
Sx that afflicts cooling how models (Ganizares, Fabbiano, & Trinchieri 1987). Gertainly 
interesting is the fate of the transient cold shell surrounding the hot bubble (especially in 
low-e solutions) during the haring activity, when the central gas surface brightness is very 
low. Due to Rayleigh-Taylor instability, the shell will break up, and perhaps cold hngers of 
gas should be observable inside the hot low density bubble, accreting on the central BH. So, 
the presence of the central heating source produces in a natural way a multiphase ISM on 
galactic scales, while the same phenomenon may be harder to obtain (Balbus 1991) in the 
cooling how scenario. 
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Fig. 5 shows the bremsstrahlung spectra in the preburst (dotted) and during burst 
(solid) phases, compared to the cooling flow (dashed) spectrum. The emitted spectrum 
is never as soft in this set of models as it is in cooling flow models, and, during bursts, 
occasionally it will have a very hard tail. 


3. Discussion and Conclusions 

In this Letter we show how the presence of a massive central BH in early type galaxies 
is able to produce naturally both the observed X-ray underluminosity with respect to the 
pure cooling flow expectations, and the large observed scatter in Lx at fixed Lb. As can 
be seen from Fig. 1, Lx - except in the very short period of bursts - is always lower than 
that of the corresponding inflow model with e = 0. Moreover, the statistical distribution 
of observed data compared with the amount of time spent at each Lx by the model here 
discussed, is eloquent (Fig. 3). Finally, due to the strong feed-back on the gas flows of 
the radiation emitted by the accretion, the total mass accumulated by the BH over 15 
Gyr is very low (~ 3 10 ®Mq), to be compared with the ~ 10 ^'^Mq of the correspondent 
e = 0-model. The same model in pure cooling-flow (without the initial SNIa driven 
wind-outflow phases) would have accumulated in its center ~ IO^^Mq of gas. 

From an observational point of view, it is interesting to note that during the accretion 
phases the galaxy luminosity is dominated by Lbh (with highest values at 10^® — erg 
s“^), while during the quiescent BH phases the galaxy emission is due only to the diffuse 
hot gas Lx. The total energy emitted by the accretion when Lbh > 10^^ is ~ 7.5 10®^ erg, 
while during the same phases the total energy emitted by the coronal gas is ~ 5.6 10®^ erg. 
The ratio between the total time spent by the galaxy when Lbh > Lx, and the total time 
spanned by the simulation is ~ 10“^: very few galaxies should be caught in a AGN-like 
phase even though most contain central BHs. 



- 8 - 

Thus, the Compton heating instability could be an alternative possibility to that 
advocated by Fabian & Rees (1995) in order to explain why the nuclei of elliptical galaxies 
are not luminous sources of radiation as expected if they host a massive central BH. A clear 
prediction of this model is that some signihcant fraction of QSOs should be embedded in 
high temperature, low surface brightness X-ray bremsstrahlung halos. 

In other works to be reported in Paper II we varied the efficiency in the range 
10“^ < e < 10“^ and the galaxy luminosity in the range 10^° < L-q/Lq < 10^^, with results 
very similar to those shown in Figs. 2-5. We are well aware that in the real accretion 
phenomenon a disk geometry for the infalling gas seems to be inescapable. Then the 
accretion luminosity will be emitted preferentially along polar directions. It is clear that (at 
least) fully 2D hydrodynamical simulations are required for a better understanding of this 
problem, and to follow the development and the hnal fate of the gas instabilities. We are 
working in this direction, but expect that many of the quantitative features of the present 
work will be carried over to the more complicated calculations. 

We would like to thank Giuseppe Bertin, James Binney, Annibale D’Ercole, Bruce 
Draine, Silvia Pellegrini, and Alvio Renzini for useful discussions and advice. This work 
was supported by NSF grant AST9424416 and in addition L.C. was supported by NSF 
grant AST9108103 and by the Italian Space Agency (ASI) with contracts ASI-94-RS-96 
and ASI-95-RS-152. 
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Fig. 1.— Panel (a): the time evolution of Lbh (bolometric) emitted at the galaxy centre. 
Panel (b): the time evolution of Lx for the model with e = 0.1 (solid line), and that of 
the same model with e = 0 (cooling flow - dashed line). Lx is calculated inside the galaxy 
truncation radius and in the range 0.5-4.5 KeV. Time interval in horizontal error bar is 
expanded in Fig. 2. 

Fig. 2.— Time expansion of the hrst burst shown in Fig. 1. The solid line is Lx, the 
dotted line Lbh- The temporal sub-structure of the burst is apparent, and the quasar-like 
luminosity 10 ^^ < LBH/ergs“^ < 10 ^^ is seen during bursts. Arrows mark epochs of “before” 
and “during” bursts referred to in text. 

Fig. 3.— The statistical distribution of Lx for observed galaxies (solid) in the range 
10.4 < log(LB/L 0 ) < 10.8 derived from Fig. 1 of CDPR. The dashed histogram represents 
the time distribution of Lx for the presented model from 9 Gyr to 15 Gyr, while the dotted 
histogram shows the cooling flow (e = 0 ) model; clearly the bursting model provides a better 
£t to the observed distribution of Lx- 

Fig. 4.— The surface brightness Sx profile (in arbitrary units) at t ~ 9.36 (dotted line, 
immediately before a burst), and t = 9.47 (solid line, during bursts), corresponding to the 
vertical arrows in Fig. 2. The dashed line is the cooling flow model at t = 15 Gyr. Note the 
very centrally bright, cuspy profile of the cooling flow, the more normal core-halo structure 
of the pre-burst profile, and the low Sx of the post-burst profile. 

Fig. 5.— Panel (a): the spectral energy distribution (in arbitrary units) of the coronal Lx, 
for the same models shown in Fig. 4. Panel (b): blow-up of the previous panel at low energy. 
Note the soft (UV) component in the cooling flow spectrum and the hard, high energy tail 
in the post burst spectrum. 
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